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Contribution to Unresolved Infrared Fluctuations from Dwarf 

Galaxies at Redshifts of 2 — 3 

Ranga-Ram Chary^, Asantha Cooray^, Ian Sullivan^ 
ABSTRACT 

In order to understand the origin of clustered anisotropies detected in Spitzer 
images between 3.6 and 8 fim, we stack the Spitzer IKAC / GieaX Observatories 
Origins Deep Survey (GOODS) images at pixel locations corresponding to faint, 
zab ~ 27 mag, optical sources with no obvious IR counterparts. We obtain a 
strong detection of the sources with a stacked median flux at 3.6/im of 130±5 nJy 
above the background. The wealth of multi-wavelength imaging data in GOODS 
enables a similar stacking analysis to be undertaken at various wavelengths be- 
tween the ultraviolet and near-infrared bands. We obtain strong stacked de- 
tections of these optically faint sources over the entire wavelength range which 
places constraints on the average properties of these sources. We find that the 
flux spectrum of the median, stacked source is consistent with a L < 0.03 L^k^uv 
galaxy with a 90% confidence interval for the redshift of 1.9 — 2.7. These sources 
produce a 3.6 /xm absolute background intensity between 0.1 and 0.35 nW m~^ 
sr~^ and the clustered IR light could account for ~ 30 — 50% of fluctuation power 
in the IR background at 4 arcminute angular scales. Although the exact redshift 
distribution of these sources is unknown, these galaxies appear to contain 5 — 20% 
of the CO- moving stellar mass density at z ~ 2.5. 

Subject headings: large scale structure of universe — diffuse radiation — infrared: 
galaxies 
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1. Introduction 

The intensity of the cosmic near-infrared background (IRB) is a measure of the total 
hght emitted by stars and galaxies in the universe which is not thermally reprocessed by 
dust. The absolute background has been estimated with the Diffuse Infrared Background 
Experiment (DIRBE; Hauser & Dwek 2001) and the Infra-Red Telescope in Space (IRTS; 
Matsumoto et al. 2005) which results in values that are about a factor of two higher than 
the intensity obtained by integrating the light from individually detected galaxies. This is 
most likely attributable to the large uncertainties associated with the removal of foreground 
zodiacal light from the DIRBE observations (Dwek et al. 2005). This hypothesis has been 
strengthened by the TeV spectrum of gamma-ray blazars which indicate that the total IRB 
intensity is smaller than the DIRBE/IRTS estimates (Aharonian et al. 2005) and that ~ 
90% of the IRB light probably arises in known galaxy populations. Yet, attempts have 
been made to explain the difference between the measured and resolved IRB intensity with 
sources at the epoch of reionization (e.g., Santos, Bromm & Kamionkowski 2002; Salvaterra 
& Ferrara 2003; Cooray & Yoshida 2004; Fernandez & Komatsu 2006). If first-hght galaxies 
are to explain the missing intensity completely, then an extreme scenario is needed with the 
conversion of at least 5% of all baryons to stars (Madau & Silk 2004). 

Instead of the absolute intensity, recent works have concentrated on spatial fluctuations 
of the IRB (Cooray et al. 2004; Kashlinsky et al. 2004). While an interpretation of unre- 
solved fluctuations is subject to extremely uncertain astrophysical modeling of underlying 
faint populations, fluctuations in deep Spitzer images have been fully attributed to first-light 
galaxies containing Pop III stars during reionization (Kashlinsky et al. 2005). An alternative 
study shows that a reasonable fraction (> 50%) of unresolved fiuctuations is arising from 
faint, unresolved sources at lower redshifts (Cooray et al. 2007). Thus, while two inde- 
pendent studies find a similar clustering amplitude for IRB fiuctuations (Kashhnsky et al. 
2004; Cooray et al. 2007), they differ in the interpretation related to the contribution from 
Pop III stars due to differences in the background hght ascribed to faint foreground galaxies 
unresolved by Spitzer. 

The suggestion for a low redshift origin in Cooray et al. (2007) comes from predictions 
using a halo model (Cooray & Sheth 2002) for the 3.6 iim population matched to their 
luminosity functions (Babbedge et al. 2006) and clustering power spectra (SuUivan et al. 
2007). In Cooray et al. (2007) this low-redshift population was identified as the faint optical 

galaxies that are resolved in the Hubble Advanced Camera for Surveys (ACS; Giavalisco et 
al. 2004) images of the Great Observatories Origins Deep Survey (GOODS; Dickinson et 
al. 2003) fields, but not detected in deep Spitzer Infrared Array Camera (IRAC) images 
of the same fields. With models, it was estimated that this faint galaxy population has a 
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3.6 fim intensity between 0.1 and 0.8 nW sr~^. The analysis in Cooray et al. (2007) 
did not exclude a high redshift contribution to the IRB, with an upper limit on the 3.6 /im 
absolute background intensity from z > 8 sources of 0.6 nW m~^ sr~^. In comparison, if 
all unresolved IR fluctuations are generated by faint z > 6.5 sources which are below the 
detection threshold of Spitzer, then the IRB intensity for such sources is >1 nW m~^ sr~^ 
at 3.6 /im (Kashlinsky et al. 2007c). 

In this Letter, we further study the physical origin of fluctuations in the IR background. 
Since faint ACS optical sources are undetected in deep Spitzer IRAC images, we stack their 
pixel locations in IRAC to establish the average IR flux. We also utilize the wealth of 
multiwavelength imaging data in the GOODS fields to obtain a stacked broadband spectral 
energy distribution of these faint optical galaxies between UV and near-IR bands. Leaving 
the redshift as a free parameter, we fit the average optical to IR flux spectrum with galaxy 
population synthesis models and find that the stacked flux spectrum is best fit by a ~ 
0.03 UV galaxy at 1.9 < z < 2.7, similar to a scaled-down population of Lyman-break 
galaxies. We also measure the expected IRB fluctuation clustering spectrum from these 
sources. 

We summarize the stacking analysis in the next section and discuss results on the average 
optical to IR flux spectrum in § 3. We discuss clustering of these faint optical sources in § 4. 
In § 5 we discuss the implications of our measurements. The magnitudes quoted throughout 
this paper are AB magnitudes. 

2. Stacking Analysis 

To establish the average IR intensity of faint optical sources that are unresolved in 
Spitzer IRAC images we stack the IRAC images at the spatial coordinates of faint, optical 
sources which are detected in the GOODS optical data. The stacking begins by first mask- 
ing all > 3(7 detected IRAC sources with mAB(3.6/im)< 26.7 mag in the GOODS/IRAC 
mosaic^. This limit is fainter than the 50% completeness limit of 24.7 mag. The mask has 
a radius of 13.5" for a source with 18 mag and scales linearly with magnitude down to 2.4" 
for sources fainter than 22 mag. We then identified sources in the GOODS ACS catalogs 
(Giavalisco et al. 2004) which are in unmasked regions of the IRAC mosaics and are there- 
fore IRAC undetected. While there are close to 22, 000 ACS sources which are unmatched 
with IRAC sources to within ~ 0.5", only 6,160 (GOODS-N) and 5,441 (GOODS-S) sources 
are sufficiently separated from brighter IRAC sources so as to be in unmasked regions of the 



description of the GOODS mosaics can be found at jhttp: / / ssc.spitzer.caltech.edu/legacy/ 
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Table 1. Results from Stacking Analysis 





GOODS-N 


GOODS-S 


Afl / arcmin^ 


169 


171 




6160 


5441 


u 


28.3 ± 0.1 


28.4 ± 0.1 


B 


28.0 ± 0.1 


27.8 ± 0.1 


V 


27.5 ± 0.1 


27.4 ± 0.1 


i 


27.2 ± 0.1 


27.1 ± 0.1 


z 


27.1 ± 0.1 


27.0 ± 0.1 


HK' 


27.3 ± 0.3 




J 




26.9 ± 0.15 


H 




26.7 ± 0.15 


Ks 




26.6 ± 0.15 


3.6 iim 


26.1 ± 0.1 


26.0 ± 0.1 


4.5 /im 


26.2 ± 0.1 


26.2 ± 0.1 


5.8 /im 


26.2 ± 0.1 


26.0 ± 0.1 


8.0 urn 


26.3 ± 0.15 


26.0 ± 0.15 



Note. — ^Acs is the total number of faint ACS sources over an effective area AD, that 
were included in the stack. Uncertainties in the stacked flux (tabulated in AB magnitudes) 
are mostly dominated by calibration systematics. 
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Fig. 1. — The average-stacked 3.6 /xm image from GOODS-N with pixel locations of all 
6160 ACS sources unmatched and unmasked in IRAC. The median stacked flux is (130 ± 5) 
nJy (or = 26.1 ±0.1 magnitude) and is detected with a statistical signal-to-noise ratio higher 
than 25. 

IRAC image. The IRAC images for each GOODS field were stacked at the pixel locations 
corresponding to these ~6000 ACS sources. The stacking was undertaken by making 12" 
image cutouts from the masked IRAC mosaic, centered on each ACS source position. The 
cutouts were then coadded using a 3a clipped median resulting in a median stacked image 
of the sources, as shown in Figure 1. A count of the number of sources which contribute 
to the stacked flux is also maintained resulting in an effective exposure map. Pixels which 
have been masked do not contribute either flux or exposure in the stacking process. Circular 
aperture photometry was performed on this stacked image and the flux measured in a 3.6" 
radius beam. The sky background was estimated from the stacked image by measuring the 
sigma-clipped median of pixels within a surrounding sky annulus spanning 5 — 6". The proflle 
of the source had a full-width at half- maximum (FWHM) of 1.6" and was consistent with 
an unresolved point source at the resolution of Spitzer. Aperture corrections were estimated 
by performing identical photometry on a calibration point source source and correspond to 
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z=2.1, M=1.4E8 Ms, A,= 1.1 mog 
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Fig. 2. — The stacked ultraviolet to IR flux spectral energy distribution of faint galaxies 
detected by ACS but not detected individually by Spitzer /IUAC in GOODS-S (black line and 
data points). The SED is consistent with the sub-L*^uv galaxy population at 2; ~ 2.5. The 
parameters of the stacked galaxy for GOODS-S are as shown in the legend. The convolution 
of the fit SED with the different passbands is shown with empty circles. The stacking in 
GOODS-N results in a best fit redshift of 2; ~ 1.9 and the SED shown as the dotted fine 
(See text for details). 

a factor of 1.14. After background subtraction and aperture corrections to the photometry 
were performed, the resultant stacked 3.6//m flux density is 130±1.3 nJy. 

In order to assess the reliability of the stacked flux, we stacked the IRAC images at 
an identical number of positions which are offset from the nominal ACS source positions by 
random amounts, up to 9" from the original source coordinates. This has the advantage of 
accounting for any possible contamination of the stacked flux due to the wings of nearby 
sources which may not be completely masked, while measuring the same local sky back- 
ground. The random stacks were repeated 100 times and the flux of the stack measured. 
The standard deviation of these values results in a flux density of 5 nJy which is a factor of 
4 larger than the uncertainty simply due to the background noise term due to source confu- 
sion and the low-level flux from the wings of masked sources. Thus, the stacking results in 
a strong detection of these faint optical galaxies and we adopt a stacked 3.6/im flux density 
of 130±5 nJy. 

This procedure was repeated at each IRAC wavelength independently for each GOODS 
fleld to assess the reliability of the result. Table 1 summarizes stacking the results for 
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all IRAC wavelengths between 3.6 fim and 8.0 /im. Through the stacking analysis, we 
obtain a strong detection of these faint optical galaxies and measure a flat flux spectrum 
(Fj, oc z/O-osiO-S), This is consistent with the flat frequency spectrum of IR fluctuations 
between 3.6 /im and 8.0 /im found by Kashlinsky et al. (2005). 

In order to estimate the total contribution to the IRB intensity of these sources, we 
multiply the median stacked flux by the exposure map. Based on the surface density of 
about 6000 galaxies in each GOODS field, we establish the absolute 3.6 /im IR intensity 
of these faint optical sources to be 0.12±0.01 nW m~^ sr~^. We note that this estimate is 
larger than the value obtained by simply multiplying the quoted number of sources (~6000) 
with the median flux. This is because the number of sources measurement requires that a 
source be unmasked out to a radius of 3.6" while the exposure map takes into account the 
exact number of sources which contribute to the stacked flux in each pixel. Since some of 
the source may be partially masked within the 3.6" radius, the number of sources is smaller 
than the number of sources as estimated from the exposure map. 

Furthermore, this estimate of 0.12 nW m~^ sr~^ only accounts for sources that remain 
outside the IRAC mask and are not confused with bright IRAC sources. Assuming the 
average flux from the stack also applies to sources affected by the mask, we estimate the 
absolute IRB intensity to be as high as 0.35 nW m~^ sr~^. This is well within the range of 
0.1 to 0.8 nW m~^ sr~^ estimated for the contribution of faint sources below the individual 
point source detection limit in IRAC images through clustering models (Sullivan et al. 2007). 

To further study the nature of these sources, beyond Spitzer/IRAC wavelengths, we 
adopted a similar procedure and stacked the same population at other wavelengths at which 
the G OODS fields have i maging data. Specifically, for GOODS-N, we utilize the KPNO U- 



band (ICapak et al. 2004f) . Hubble/ ACS BViz (Giavalisco et al. 2004), and near-infrared HK' 
JCapak et al. 2004h . For GOODS-S, we use the CTIO U-band and the ESO/ISAAC near- 
infrared data in the JHK bands in addition to the Hubble and Spitzer data. For the Hubble 
data, we simply measure the median flux values of the sources since they are individually 
detected at the ACS wavelengths. For the other wavelengths, we excised subimages for 
each source and estimated a sigma clipped average of the subimages to obtain the stacked 
flux, as was done for the IRAC data. Aperture corrections, background estimates and noise 
calculations were performed as described for the IRAC data. 
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Fig. 3. — The top solid circles show the angular power spectrum of unresolved IR fluctuations 
in IRAC channel 1 (3.6 fim) of GOODS-N field. The solid squares show the angular power 
spectrum that masks the locations of faint ACS sources that are undetected in IRAC. The 
same power spectrum published in Cooray et al. 2007 is shown with open squares and the 
difference accounts for a new correction we have introduced here for the window function of 
the ACS source mask. With ACS sources masked, the fluctuation power spectrum amplitude 
is reduced at £ ~ 7 x 10'^ by a factor of about 2 (±0.6). The circles and dashed lines below 
show the predicted power spectrum for faint optical sources with their 3.6 fim number 
counts distributed with two faint-end slopes: the open and filled circles are for clustering 
in GOODS-N and GOODS-S assuming the steeper slope for counts with a ~ 0.6, while 
the dashed lines bracket the clustering for a flatter slope with a ~ {dN/dm oc m" when 
m > 26) For reference, the two solid lines are predictions from clustering models of the faint 
source population with the range covering the uncertainty in model parameters (Sullivan et 
al. 2007). 



3. Redshift and Stellar Mass Density Estimates 

In order to understand the physical properties of these faint galaxies, we fit the stacked 
flux in Table 1 with Bruzual & Chariot (2003; BC03) population synthesis models. Redshift, 
mass, extinction, age and e-folding time of star-formation were left as free parameters. The 
redshift range was 0.1 — 6 in steps of 0.1, extinction range was Ay = — 5 in steps of 
0.1, age t, varied from 10^ yr to the age of the Universe at each redshift, while the e- 
folding time r spanned the entire range from instantaneous to constant star-formation. For 
GOODS-N, we find a best fit redshift of 1.9, stellar mass of 1.3x10*^ Mq, Av=1.4 mag, t=lO 
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Myr and r=200 Myr, indicating an ongoing reddened starburst (Figure 2). For GOODS- 
S, which has far better near-infrared data, we find a best fit redshift of 2.1, stellar mass 
of 1.4x10^ Mq, Av=1.1 mag, t=15 Myr, r=200 Myr. The 90% confidence interval for 
the redshifts as determined from the lowest chi-square fit at each redshift, spans the range 
1.9 < z < 2.9 for GOODS-S and 1.3 < z < 2.7 for GOODS-N. The [/-band detection 
implies that the majority of the sources must be at 2; < 2.5. The sharp increase in the flux 
in the IRAC channels compared to the NIR bands implies that the 1.6/im bump must be 
in the IRAC passbands. This implies z > 1.3. Thus, even a qualitative comparison of the 
stacked flux densities with a typical galaxy SED results in a redshift range similar to what 
we estimate from the chi-square values. We adopt 1.9 < z < 2.7 as the redshift range of 
these sources. The uncertainty in the derived stellar mass is about a factor of two, ranging 
from 1.3 - 3.0 x 10^ Mq. 

For our best fit template, the rest-frame 1500A UV luminosity of the stacked galaxy is 
1.8x10^ L0 while its rest-f rame V-band lumin o sity is 1.2 x 10^ Lg). Sinc e L=k^uv,z=3=5.8x10^° L© 



and L^,^v,z=3=8xl0^° Lq (ISteidel et al. 19991 : iMarchesini et al. 20071 ) . it implies that these 



faint objects are galaxies which are about ~ 30 — 60 times fainter than the characteristic 
luminosity of field galaxies at these redshifts. 

Although the exact redshift distribution of these faint galaxies is not known at this time, 
it is illustrative to estimate the fraction of the stellar mass density hidden in these objects. 
We assume that the ~6000 galaxies that we have stacked on are distributed between 1.9 < 
z < 2.7 with a stellar mass of 1.4x10® Mq. This corresponds to a co-moving stellar mass 
density of 2x10^ M© Mpc'^ which is a strong lower limit. Assuming that our stacked flux is 
typical for the ~22000 ACS sources which are unmatched to IRAC sources, would imply a 
total stellar mass density in faint galaxies of 6.7x10^ Mq Mpc~'^ at 2; ~ 2.4. For comparison, 
the stellar mass densit y in Lyman-break galaxies at these redshifts is 3.6x10'' Mq Mpc~^ 



(jPickinson et al. 20031 ). Thus, the faint, sub-L* galaxy population could account for upto 



20% of the total stellar mass density at z ^ 2.4. 



4. Expected Clustering 

We also measure the clustering of these sources to establish the amplitude of unresolved 
fluctuations at IRAC bands. Since we are interested in the anisotropy of IR light, to establish 
the power spectrum of fluctuations, we need the flux distribution within this population. The 
stacking analysis only allows us to establish either the total or the average flux of this sample. 
This forces us to make an estimate of the flux distribution and use that to predict clustering 
of unresolved fluctuations produced by these sources. The angular power spectrum of IR 
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fluxes associated with these sources is determined using the same technique as described in 
SuUivan et al. (2007). We account for the flux distribution by randomly assigning IR flux 
to the ~ 6000 optical source locations between 70 nJy and 700 nJy such that the total flux 
assigned is the same as the total flux measured for these sources from the stack. We assign 
fluxes such that the number counts trace the extrapolated slope from known IRAC counts 
down to mAB(3.6;um) of 24.7 mag. Given the uncertainty in the faint-end slope, we take 
the slope a. [dN/dm oc rrf) to be either 0.6 or 0.0. We tested how our predictions change 
with variations to these parameters and found consistent results within errors. Since the 
flux assignment is random, to get an average of the expected clustering, we randomize the 
assignments and the fluxes and use a Monte-Carlo approach to obtain the mean and variance 
of clustering. 

We summarize our results in Figure 3, where we also compare with a direct measurement 
of the clustering of IR fluctuations and expectations based on the clustering models of Sulli- 
van et al. (2007; solid lines). In Figure 3, we also show the IR fluctuation spectrum when the 
faint optical sources that we stack on are masked from the IRAC image with open squares 
(Cooray et al. 2007). The solid squares show the revised measurement after correcting for a 
bias associated with the window function introduced by the large mask. To understand this 
correction noted in Kashlinsky et al. (2007a), we note that the fluctuation spectrum measur- 
able in Fourier (multipole) space in the presence of a mask is Cg, = J2e' ^U'Ce, where Ww 
is the window function associated with the mask and Q/ is the power spectrum of interest. 
We recover the latter with measurements of Q by flrst generating associated with our 
mask (Appendix A of Hivon et al. 2002) and then iteratively inverting for Q using an inver- 
sion similar to Dodelson & Gaztaiiaga (2000). The inversion agrees to 10% with a separate 
estimate of the angular power spectrum using a likelihood method where one maximizes the 
likelihood of an estimated C^^* with a model spectrum C™^ using C^^* = W^^/C™^ given 
the measurements and the error a^. This procedure can be described as minimizing the 
with x' = {Ce - Crf/al 

As shown in Figure 3, the difference between the fluctuation amplitude of unresolved IR 
fluctuations and the fluctuation amplitude with ACS sources masked is about a factor of 2 
(± 0.6) at £ ~ 7 X 10^. Our measurements then suggest that, from this difference, up to 50% 
of IR fluctuations can be accounted by faint ACS sources. It could be that by accounting 
for further fainter optically sources, this fraction can be increased, but it is quite unlikely 
that 100% of IR fluctuations reported in Kashhnsky et al. (2005, 2007b) are generated by 
z > 6.5 galaxies hosting mostly Pop III stars that are optically invisible. 

The difference between IR fluctuations with and without ACS sources masked must 
be reproduced by the clustering of IR hght produced by the same ACS sources. With 
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flux assignments based on the stacking analysis, we find that the difference is reproduced 
(though these indirect estimates of the expected clustering IR fluctuations produced by faint 
ACS sources are uncertain due to uncertainties in the flux distribution). We conclude that 
irrespective of whether we measure the power spectrum in the IRB fluctuations after masking 
optical sources or whether we measure the power spectrum of faint optical galaxies with an 
average flux based on our stacking analysis, we find that the power in faint {zab ~ 27 
mag), z ~ 2.5, galaxies accounts for at least ~30— 50% of the power in the IR background 
fluctuations on angular scales of ~4', where the measurement is not strongly sensitive to 
cosmic variance. 

On similar angular scales, Kashlinsky et al. (2007a) find a power of 0.015 nW m~^ sr~^ 
for the optical galaxies with zab >26.6 mag in HDFN-E2, without a quantification of the un- 
certainty. They also demonstrate that the power increases as brighter galaxies are included. 
40% of the galaxies that go into our stacking analysis are brighter than this threshold and 
therefore the power they measure is a lower limit to our power spectrum shown in Figure 3. 
This can be compared to their measured value of 0.035 nW m~^ sr~^ for the amplitude of 
power in the IR background fiuctuations in the HDFN on the same angular scale/multipole 
moment. Making use of their estimate of the fiuctuations associated with faint, optical 
galaxies at low redshifts, we find that Kashlinsky et al. (2007a, 2007b) also find a ~ 40% 
contribution from such sources to fluctuations at ~ 4 ' angular scales, though their contri- 
bution is ignored in their subsequent interpretation that assumes all fluctuations are from 
Pop III stars before reionization. The analysis of the optical sources in GOODS-S is not 
presented in Kashhnsky et al. (2007a). The fundamental discrepancy between our two 
groups arises on larger angular scales, at 400 arcsec or larger, where they claim the power 
spectrum is rising. We are unable to constrain the power spectrum on such large angular 
scales {£ ~ 1500) due to the cosmic variance associated with the limited size of the GOODS 
fields used in our analysis, which becomes important at angular scales above 400 arcseconds. 



5. Discussion 

The Kashlinsky et al. (2007c) interpretation for IR fiuctuations involves z > 6.5 galaxies 
with an IR background intensity of > 1 nW m~^ sr~^ at 3.6 fim. This is a large intensity 
given that all resolved sources, so far, lead to an IR background intensity between ~ 6 nW 
m~^ sr~^ (Sullivan et al. 2007; Fazio et al. 2004) and ~ 10 nW m~^ sr~^ (Levenson & 
Wright 2008). To avoid individual detections of a large number ol z > 6.5 galaxies, such a 
high background intensity must then be hidden in a source population with a large surface 
density but with individual faint fluxes of around 10 nJy for each source. Given the large 
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surface density of the population, pixel to pixel intensity variations are a small fraction of 
the total intensity produced by those sources. 

The alternative explanation is that a reasonable fraction of fluctuations arises from a 
population of sources that have a low surface density but are just below the <5pzteer/IRAC 
detection threshold. The results from the stacking analysis presented here suggest that 
faint optical galaxies, with an average source flux of around 130 nJy at 3.6 //m, contribute 
between ~30— 50% of the power in the fluctuations on the largest angular scales that we are 
able to measure. The average IR absolute intensity produced by these sources is in the range 
0.12-0.35 nW m~^ sr~^ at 3.6 /im, in agreement with the Kashlinsky et al. (2005) estimate 
that the contribution from faint extragalactic sources is ~0.15 nW m^^ sr^^. However, they 
claim that such faint objects cannot account for the strong clustering signal in unresolved 
IR light. We have demonstrated in this paper that a significant fraction of the fluctuations 
do arise from these faint galaxies. These sub-L* galaxies contribute up to 20% of the stellar 
mass density at 2; ~ 2.5 and can be described as a low luminosity version of well-studied 
Ly man-break galaxies. 

While these sources do not fully explain all IR fluctuations measured with deep IRAC 
images, in Cooray et al. (2007), we placed a conservative upper limit that any contribution 
from z > 6.5 sources must have an absolute intensity less than about 0.6 nW m~^ sr~^. 
Beyond Spitzer IRAC, studies have also been conducted at lower near-IR wavelengths with 
NICMOS. Though limited by the small fleld-of-view of NICMOS, they show that IR fluc- 
tuations at 1.6 and 1.25 /xm are more consistent with a 2; < 8 origin than a high redshift 
interpretation (Thompson et al. 2007a, 2007b). Beyond IR fluctuations, it will be inter- 
esting to further study the nature of the faint, dwarf galaxy population at 2; = 2 to 3 to 
understand the contribution it makes towards bridging the systematic offset between the 
co-moving stellar mass density and the integrated star-formation rate density (Hopkins & 
Beacom 2006). 

Acknowledgments: We thank Mark Dickinson for helpful suggestions. We also acknowl- 
edge the contributions of various members of the GOODS team who are associated with 
key aspects of data processing and handhng. This research was funded by NASA grant 
NNX07AG43G, NSF CAREER grant AST-0645427, award number 1310310 from Spitzer 
for Archival Research, and program number HST-AR-11241.01-A by NASA through a grant 
from the Space Telescope Science Institute, which is operated by the Association of Univer- 
sities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555. 



- 13 - 



REFERENCES 

Aharonian, F. et al. 2006, Nature, 440, 1018 

Babbedge, T. S. R., Rowan-Robinson, M., Vaccari, M. et al. 2006, MNRAS, 370, 1159 

Bruzual, G., & Chariot, S., 2003, MNRAS, 344, 1000 

Capak, R, et al., 2004, AJ, 127, 180 

Cooray, A. & Sheth, R. 2002, Phys. Rep., 372, 1 

Cooray, A. et al. 2004, ApJ, 606, 611; erratum-ibid. 2005, ApJ, 622, 1363. 
Cooray, A. & Yosliida, N. 2004, MNRAS, 351, L71 
Cooray, A. et al. 2007, ApJ, 659, L91 

Dickinson, M. et al. 2003, in Proceedings of the ESO Workshop, The Mass of Galaxies at 
Low and High Redshifts", eds. R. Bender & A. Renzini, p. 324 ( |astro-ph/ 0204213) 

Dickinson, M. et al. 2003, ApJ, 587, 25 

Dodelson, S. & Gaztaiiaga, E. 2000, MNRAS, 312, 774 

Dwek, E., Arendt, R. & Krennrich, F. 2005, ApJ, 635, 784 

Giavalisco, M. et al. 2004, ApJL, 600, 93 

Fazio, C. G. et al. 2004, ApJS, 154, 39 

Fernandez, E. & Komatsu, E. 2006, ApJ, 646, 703 

Hauser, M.G. & Dwek, E. 2001, ARA&A, 39, 249 

Hopkins, A. M., & Beacom, J. F., ApJ, 651, 142 

Hivon, E. et al. 2002, ApJ, 567, 2 

Kashlinsky, A. et al. 2004, ApJ, 608, 1 

Kashlinsky, A. et al. 2005, Nature, 438, 45 

Kashlinsky, A. et al. 2007a, ApJ, 666, LI 

Kashlinsky, A. et al. 2007b, ApJ, 654, L5 



-14- 



Kashlinsky, A. et al. 2007c, ApJ, 654, LI 

Levenson, L. & Wright, E. L. ApJ submitted, larXiv:0802.T239l 

Madau, P. & Silk, J. 2005, MNRAS, 359, L37 

Marchesini, D., et al. 2007, ApJ, 656, 42 

Matsumoto, T. et al. 2005, ApJ, 626, 31 

Salvaterra, R. & Ferrara, A. 2003, MNRAS, 339, 973 

Santos, M. R., Bromm, V., & Kamionkowski, M. 2002, MNRAS, 336, 1082 

Steidel, C. C, et al. 1999, ApJ, 519, 1 

Sullivan, I. et al. 2007, ApJ, 657, 37 

Thompson, R. et al. 2007a, ApJ, 657, 669 

Thompson, R. et al. 2007b, ApJ, 666, 658 

Wright, E. L. & Reese, E. D. 2000, ApJ, 545, 43 



This preprint was prepared with the AAS I^TJrjX macros v5.2. 



